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The thermoelectric properties of La-filled skutterudites are discussed from the point of view 
of their electronic structures. These are calculated from first principles within the local density 
approximation. The electronic structure is in turn used to determine transport related quantities. 
Virtual crystal calculations for La(Fe,Co)4Sbi2 show that the system obeys near rigid band behavior 
with varying Co concentration, and has a substantial band gap at a position corresponding to the 
composition LaFe3CoSbi2. The valence band maximum occurs at the V point and is due to a singly 
degenerate dispersive band, which by itself would not be favorable for high thermopower. However, 
very flat transition metal derived bands occur in close proximity and become active as the doping 
level is increased, giving a non-trivial dependence of the properties on carrier concentration and 
explaining the favorable thermoelectric properties. 



There has been a recent revival of activity in the search 
for improved thermoelectric (TE) materials, with an em- 
phasis on novel materials concepts!! The utility of a ma- 
terial for TE applications is measured by its dimension- 
less figure of merit, ZT = (tr / k)S 2 T ', where T is temper- 
ature, a is the electrical conductivity, S is the Seebeck 
coefficient and k is the thermal conductivity, which con- 
tains both electronic and lattice contributions, n e and 
Ki, respectively. High values of S are characteristic of 
highly doped semiconductors; unfortunately, these have 
low electric conductivity and thus small ratio (a/n). On 
the other hand, metals usually have this ratio close to 
the Wiedemann- Franz value, a/n e . If the Wiedemann- 
Franz law holds, S — 160 V/K is required for the state- 
of-the-art value ZT = 1. Much TE research over the 
past 40 years has focussed on covalent semiconducting 
compounds and alloys, composed of 4th and 5th row el- 
ements, with a view to finding low thermal conductiv- 
ity materials that have reasonable carrier mobilities and 
high band masses. Most current generation TE materi- 
als are of this type, e.g. Bi 2 Te3, Si-Ge, PbTe. Despite 
research efforts spanning three decades, little progress 
in increasing ZT has been achieved, and in particular 
Bi2Te3/Sb2Te3 has remained the material of choice for 
room temperature applications. 

Recently three new materials with ZT > 1 have been 
reportecoQ, and these do not clearly fall in the same 
class as previous TE materials. fi-Z^Sbz, with reported 
ZT up to 1.3 has a large region of linear temperature de- 
pendence of the resistivity, suggestive of a metallic rather 
than semiconducting material, but unlike normal metals 
this is accompanied by high thermopowers. CeFe4Sbi2 
both by itself and doped with Co shows high values of ZT 
combined with very low thermal conductivities, and de- 
pending on the composition can show either metallic-like 
or semiconducting temperature dependencies of the re- 
sistivity. At low temperature, CeFe4Sbi2 displays prop- 
erties reminiscent ofpheavy fermion materialsd. Band 
structure calculationsQ show that the Ce /-states, indeed, 
contribute significantly to the electronic structure near 



the Fermi energy, leading the enhanced band masses, fa- 
vorable for TE, and producing band gaps via hybridiza- 
tion with the valence states. This would imply a rather 
different electronic structure for LaFe4Sbi2. However, 
La(Fe,Co)4Sbi2 does have ZT w 1 for appropriate con- 
ditions. Here we report first principles calculations for 
La(Fe,Co)4Sbi2 within the local density approximation 
(LDA), similar to our previous calculations for binary 
and Ce-filled skutteruditesaia. We then compute trans- 
port properties, based on the calculated band structure, 
and show them to be in accord with the experiments. 
An analysis of the underlying band structure reveals the 
mechanism for high thermopower, and suggests ways to 
optimize the TE properties. 

The calculations were performed using the gen- 
eral potential linearized augmented planewave (LAPW) 
methodd. This method makes no shape approximations 
to either the potential or charge density and uses flexible 
basis sets in all regions of space. As such it is well suited 
to materials with open crystal structures and low site 
symmetries like La(Fe,Co)4Sbi2. Well converged basis 
sets of «2100 functions were used. Local orbital exten- 
sions to the basis set were used in order to relax lineariza- 
tion errors generally and to include the upper core states 
of La consistently with the valence states. 

The calculations were based on the experimental crys- 
tal structure of LaFe4Sbi2. The electronic structures of 
La(Fe,Co)4Sbi2 alloys were calculated as for LaFe4Sbi2 
using the virtual crystal method, i.e. self-consistent 
LAPW calculations were performed using average ions 
with fractional charges intermediate between Fe and Co. 
This approximation includes the average ionic charge and 
band filling, but is known to distort the electronic struc- 
ture for alloys when the scattering properties of the ions 
at the mixed site differ substantially. Calculations were 
performed for 0%, 25% and 50% Co substitution. The 
calculated bands near the Fermi level differ little, provid- 
ing aposteriori justification of the virtual crystal approxi- 
mation. Transport properties relevant to TE were deter- 
mined from the calculated band structures using stan- 
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dard kinetic theory as given by Ziman and other, 
a(T) = ^JdeN(e)v 2 (e)r(e,T) (-^M 

1 f ' ( T s ( df(e) 
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Here N(e) is the density of electronic states at the energy 
e per unit volume, r is the scattering rate for electrons, 
and the average velocity v(e) of electrons with the energy 
e is defined via 



N{e)v\e) = {2tt)- a / <J(e k - e)dk 
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The above formulae are readily calculated from the 
band structure, provided that the scattering rate r is 
known. Fortunately, in most cases r does not vary too 
much with the energy, although there are exceptions to 
this rule, e.g. Pd metalO where Ep occurs near a very 
sharp feature in N(e\- and Kondo systems where there 
is resonant scattering^. La(Fe,Co)4Sbi2, on the other 
hand, shows no indications of such behavior, and r(e) 
can be replaced by a constant which then cancels out 
in many transport properties, like Seebeck coefficient or 
Hall number. 




FIG. 1. Virtual crystal band 

structure of La(Feo.75Coo.2s)3Sbi2 in the vicinity of the 
Fermi level. 

The calculated band structure of La(Feo.7sCoo.25)3Sbi2 
near the Fermi level is shown in Fig. |l|. Although t 
initial idea of alloying unfilled skutterudites with L 
was to improve the figure of merit by reducing the lattice 
thermal conductivity, it turns out that La plays a posi- 
tive role from the point of view of the electronic structure 
as well. First principles calculations for CoSb30 show a 



small band gap with a highly non-parabolic, quasi-linear 
valence band dispersion.—. This result, which now has 
experimental supportE3lla, combined with the kinetic 
transport theory above, implies that in p-type CoSb3 
the doping level, n, dependence of S and a will differ 
from that of a semiconductor with parabolic bands. In 
particular, in the constant scattering time, degenerate 
regime, S oc n" 1 / 3 and a oc n 2 / 3 instead of the usual 
S oc n 2 / 3 and a oc n. This means that the power fac- 
tor, S 2 will be less dependent on doping level than in 
a conventional case and that it will be more difficult to 
optimize ZT by adjusting n (but note that k has an 
electronic component proportional to a). 

The calculated band structure of La(Feo.75Coo.25)3Sbi2 
shows an indirect T — N gap of 0.60. This is qualitatively 
different from the band structure of CoSb3. The direct 
gap at r is 0.76 eV. Both the conduction and valence 
band edges are formed from parabolic bands with mostly 
Sb p— character. In addition there are some rather flat 
(heavy mass) primarily Fe/Co d— derived bands near but 
not at the band edge both above and below the gap. 

To understand this band structure we first look at the 
crystal structure. A useful starting point is a hypothet- 
ical LaFe compound with the CsCl (B2) structure. We 
then construct 2x2x2 supercell, keep two La atoms along 
(111) directions, and remove the remaining six La's. The 
resulting pores are filled with the Sb4 rings, oriented in 
such a way as to keep the (111) threefold axis. The Sb4 
rings form quasimolecules, since the Sb-Sb bonds are the 
shortest in the crystal (« 2/3 of the Fe-Fe bond length). 
12 p— states of Sb atoms in a molecule hybridize with each 
other forming four strongly bound ppa bonding states 
(split into two close pairs by weak ppn interaction), four 
ppa antibonding states, and two pairs of weakly bound 
ppn bonding/antibonding states, formed by the orbitals 
perpendicular to the plane of the ring. Half of these states 
are occupied by 6 p— electrons of four Sb atoms. Fe/Co 
d— states are close in energy to the Sb p-states. They 
give rise to two sets of bands derived from 6 ti g states at 
lower energy and 4 e g states at higher energy per metal 
ion. Three Fe and one Co together have 33 d— electrons, 
of which 24 reside in the t^g bands and 9, together with 
three La electrons, fill up twelve more states in Sb rings. 
Thus, the highest occupied states in the ionic picture oc- 
cur between antibonding ppix Sb bands and ti g Fe/Co 
states. Let us first discuss the former. The symmetry 
of the relevant combination of p— orbitals with respect to 
the center of the ring is xyz, that is, an /—symmetry. 
Not surprisingly, these combinations hybridize with the 
La xyz states, forming a quasi-//<5 band with the max- 
imum at the r point, where the ffS interaction is the 
strongest (by symmetry). The symmetry of LaFe4Sbi2 
is I3m and does not have the full set of 48 cubic opera- 
tions, but the effective mass of this specific band is very 
isotropic, and rather small, m* « 0.2m e . 

Let us now consider the second highest occupied band 
at the r point, which is a Fe/Co t^g band. With maximal 
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cubic symmetry this band would have been triply degen- 
erate at r. Since the actual symmetry of LaFe4Sbi2 is 
lower, the corresponding state is split into double de- 
generate and nongencrate states, the latter being the 
most symmetric combination of the three t^g orbitals, 
xy + yz + zx. This state can also be written as 3z' 2 — r 2 , 
where z 1 is the threefold axis, z' is one of the three prin- 
cipal axis for the mass tensor at the T point. Note that 
this state is directed towards La and thus has little over- 
lap with the d orbitals on the neighboring Fe/Co sites. 
This results in a low dispersion and high effective masslij, 
(m*) « 3. The top of this band is only 0.1 eV below the 
top of the previously described Sb band. The combina- 
tion of the two bands of so different effective masses in 
direct vicinity of the Fermi level is what provides high 
thermopowers in this material. 




E-Ep, mRy 



FIG. 2. Hole concentration (+), Hall number (dashed) 
and thermopower at 600 K (solid) as a function of the 
Fermi level shift with respect to the center of the gap 
(Fermi level for stoichiometric La(Feo. 75 Coo. 25)38012. 

In order to understand the transport more quanti- 
tatively, we have performed calculations of transport 
propertiestll as a function of hole doping away from the 
gap based on the 25% Co virtual crystal band structure of 
Fig. |. The Hall number, carrier concentration (doping 
level), and thermopower at 600K are shown as a function 
of Ep in Fig. |^. Over almost all of the range shown the 
Hall number is lower than the doping level. This devi- 
ation, which increases with band filling, reflects devia- 
tions from isotropic parabolic band behavior. Although 
the presence of the heavy band is responsible for the high 
thermopower of this material, S at 600K does not show 
any noticeable structure near —26 mRy, which is the on- 
set of the heavy band. This is because the form of Eq. (Q) 
has strong contributions at energies up to 2-3 kT from 
the chemical potentialta and because kT at 600 K is 3.7 
mRy. Thus the heavy band contributes to S over the 
entire range of Fig.^|. However, the carrier concentration 
increases sharply below the onset of the second band, re- 
flecting its flat dispersion. This provides a mechanism for 
pinning the Fermi level near the onset, providing high 



thermopowers even though the carrier concentration is 
apparently difficult to control in these materials. 

Although the Hall number is not reported^ the 
sample shown by Sales, Mandrus and Williamsu has 
5(600K)«180 V/K. Comparing with our calculations, 
this corresponds to a chemical potential of —30 mRy, or 
4 mRy below the onset of the heavy band. This implies 
a doping level of 0.08 holes per cell and a Hall number 
of 0.045 holes per cell. The corresponding Hall concen- 
tration of 1.2 x 10 20 holes/cm -3 is consistent with the 
reported lower end of the range of Hall concentrations 
reported in Ref. [| for various samples (1 x 10 20 cm-3). 




FIG. 3. Calculated at E - E F = 30.3 mRy (dashed) 
and experimental (solid) thermopower. 

The temperature dependence of S at a chemical po- 
tential —30.3 mRy (matching the experimental value of 
S at 500K) is shown in Fig. Q. The agreement is fair; the 
slope at low temperatures is too high and the saturation 
temperature is too low, although the qualitative fact that 
S is rather high even at high doping levels is reproduced. 
One possible reason for the quantitative disagreement is 
neglect of the energy dependence of r. Since the two rele- 
vant bands have very different characters (one is Sb p and 
another is Fe/Co d), even if the relaxation rates in each 
band are energy independent, as is likely, they may be 
different between the two bands, yielding an effective en- 
ergy dependence of r. This can change noticeably the 
temperature dependence. Roughly speaking, the con- 
ductivity is defined mostly by the light band, thus the 
thermopower is more or less additive for the two bands. 
Because of very different dispersion laws the tempera- 
ture dependence of partial thermopower in each band is 
quite different, thus changing the balance between the 
two can also change the net temperature dependence for 
the total thermopower. There is, however, another, per- 
haps more important effect, namely that the calculated 
band structure is distorted relative to the samples due to 
the assumption of full La filling. We note the fact that 
the actual, as measured, material is La deficient, and we 
observe from our calculations that La f xyz orbitals are di- 
rectly involved in the formation of the highest occupied 
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( "light" ) band. The relative positions of the heavy and 
light bands, which are the key ingredient in determining 
the temperature dependence of S, thus depend strongly 
on La filling. 

Related to this it should be noted that the thermal con- 
ductivity for this sample in the temperature range above 
500 K is roughly 1/3 electronic and 2/3 lattice in origin. 
As mentioned Ki is strongly reduced from ni of CoSb3 
by La addition, but most of this effect is expected to 
occur at fairly low La concentrations within the rattli&g 
ion framework of Slack and Cahill, Watson and Pohfca, 
and supported by measurements in the Ce(Ee,Co)4Sbi2 
on samples with different Ce concentrations!! Because 
of the strong interaction between the valence bands and 
La, La vacancies should strongly scatter electrons reduc- 
ing the electrical conductivity. This is consistent with the 
wide range in hole mobilities (2-30 cm 2 /V-s) measured 
for the various samples. Based on this we conjecture that 
samples with higher La filling, if they can be made, would 
have higher mobilities and higher values of a at a fixed 
band filling. Because of the Wiedemann-Franz relation, 
this would also lead to higher K e . In that case, the value 
of S would become the most important factor determin- 
ing ZT at a given temperature. We speculate it may be 
quite feasible to increase ZT significantly in this material 
if the La concentration could be made stoichiometric near 
25% Co concentration. In this case, the doping level for 
maximum ZT will occur at lower carrier concentrations 
than any of the reported samples. Depending on how 
much the mobility can be improved by La filling, it may 
be possible to obtain high values ZT at lower tempera- 
tures by doping so that the chemical potential is roughly 
kT above the onset of the heavy band (i.e. in the 10 19 
cm~ 3 range). 

To summarize, we have reported first principles cal- 
culations of band structure and transport properties 
of La(Fe,Co)4Sbi2. Despite high ZT and large ther- 
mopower this material may be reasonably described as 
metallic, both from the point of view of experimental 
measurements and the calculated band structure. The 
electronic structure of La(Fe,Co)4Sbi2 is strongly dis- 
torted from that of the corresponding binary, CoSb3 im- 
plying an important role of La in formation of the valence 
bands. This fact finds its explanation in the specific char- 
acter of the Sb states that form the highest valence band. 
Our results suggest that samples with higher La filling, 
and lower hole concentrations may have even better TE 
properties. 
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